Magnetic strong coupling in a spin-photon system and transition to classical regime 
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We study the energy level structure of the Tavis- Gumming model applied to an ensemble of 
independent magnetic spins s — 1/2 coupled to a variable number of photons. Rabi splittings 
are calculated and their distribution is analyzed as a function of photon number nmax and spin 
system size N. A sharp transition in the distribution of the Rabi frequency is found at Timax ~ A'^. 
The width of the Rabi frequency spectrum diverges as \/iV at this point. For increased number 
of photons nmax > N, the Rabi frequencies converge to a value proportional to ^nmax- This 
behavior is interpreted as analogous to the classical spin-resonance mechanism where the photon is 
treated as a classical field and one resonance peak is expected. We also present experimental data 
demonstrating cooperative, magnetic strong coupling between a spin system and photons, measured 
at room temperature. This points towards quantum computing implementation with magnetic spins, 
using cavity quantum-electrodynamics techniques. 

PACS numbers: 42.50.Ct, 71.45.-d, 75.45. +j, 64.60.-i 



I. INTRODUCTION 



Interactions of quantum systems with electromagnetic 
excitations are at the core of quantum information pro- 
cessing. Using photons and photonic entanglement, 
qubits can be detected and manipulated, and quantum 
information can, in principle, be transferred over long 
distances^. Of particular interest are the resonant modes 
in electromagnetic cavities, which have the potential of 
inducing a strong coupling regime such that the interac- 
tion outlast both photon's decay and qubit decoherence 
times^''^. Following the work of Dicke^ on multi-atom su- 
perradiance, the case of a single atom in interaction with 
n photons has been studied theoretically by Jaynes and 
Cummings^, and later on generalized®^^ for a number of 
N otherwise non-interacting atomic systems. Other the- 
oretical studies, applied to solid state systems^°, have 
included environmental effects as well (e.g., in semicon- 
ducting materials^^'^^) or ensemble-locking in a giant 
spin^'^. Experimentally, the phenomena of strong cou- 
pling regime has been reached by using the electric field 
component of the electromagnetic excitations: in one or 
more atomic systems^^'^^, semiconductors^® and super- 
conducting qubits in interaction with one^^ or more pho- 
tons"'^®'"'^^. These studies prove the appearance of the so- 
called vacuum-field Rabi splitting (VRS) in the absorb- 
tion peak of a probing photon field. 

In contrast, achieving large magnetic coupling between 
a photon and a quantum spin, has been explored to a 
lesser extent, due to the typical smallness of the mag- 
netic component (B-field) of the electro-magnetic field. 
However, since spin-based qubits do reveal significant co- 
herence times for temperatures up to ambient value^®'^^, 
the issue of coupling spin qubits to photons for data ma- 



nipulation and transfer becomes of increasing interest. 
In the usual magnetic resonance methods, e.g. electron 
spin resonance (ESR), the absorbtion measurement of the 
electromagnetic field is related to the energy structure of 
the spin system. Feedback effects of the B-field compo- 
nent on the spins are ignored which means that for a two- 
level system there is one absorbtion peak at a frequency 
matching the levels separation. On the other hand, as 
mentioned above in the case of electrical coupling, pho- 
ton absorption can probe the quantum mechanical in- 
teraction between the quantum system and the cavity 
photons which leads to VRS. In this paper, we will study 
the relation between these two cases, one classical and 
the other one quantum, by comparing the effects of mag- 
netic coupling between N non-interacting spins s — 1/2 
and the external radiation field. The transition between 
the classical and quantum case will be analyzed as well. 

Because the wave length of the external field is large 
compared with the distances between spins, all the spins 
interact with a single mode of electromagnetic field. In 
the classical case of spin resonance, the system Hamil- 
tonian is given by Hs = "Hzo + ^ac where Hzo is the 
Zeeman coupling to a static field H^'. 



N 



(1) 



where niz = —gsfJ.sO'i /2 gives the magnetic moment of 
spin i (Qs is the g-factor, equal to 2 for a free spin and 
/is is the Bohr magneton) and hujo is the Zeeman split- 
ting generated by Hz. The term Hac represents the 
spin coupling to an alternating B-field component ho 
oscillating with a frequency a;/27r. Using the notation 
hflji = fiohogsiiB /'^ and the Pauli projection and rais- 
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ing/lowering operators this term is written as: 

1 ^ 

HacR = ^n^n ^ (e-V+ + Vr) , (2) 
for a rotating B-field or, for an uniaxial B-field, as 

N 

^acx = nflR cos(a;t) ^ . (3) 

i=l 

When treating the radiation field quantum mechani- 
cally, as in the Jaynes-Cummings model, the spin-photon 
coupling is described by a parameter g, here assumed to 
be the same for all spins: 

N 

n^Hzo + hgJ2 (K^ + b^'^i) + hojb^b, (4) 

i 

with b,b^ the photon annihilation/creation operators. 

The amplitude of the electro-magnetic field is a dynami- 
cal variable, but not an external parameter as in the case 
of Eq. (2) or (3). Namely, the strength of the electro- 
magnetic field is given by b^b in the quantum mechanical 
model in Eq. (4), while it is given by ho in the classical 
models in Eqs. (2) and (3). 



Rabi oscillation. By transforming the wave function 

\m) as 

\^{t)) = 6'^"°*"' mt)) = J7|$(i)), (6) 
the rotating frame version of Hs is given by 

^ N N 

UHsU-^ = -nnn ^ (a+ + a") = ^ a„ (7) 

i i 

at resonance. In this representation, the z component 
of the magnetization rotates around the x-axis with the 
angular velocity 9 = Qr, which is the Rabi oscillation. It 
should be noted that the phase oscillates with the angu- 
lar frequency cjq i which causes a rotation of the magne- 
tization around the z axis in the laboratory frame. The 
effective eigenvalues of Hamiltonian (7) are given by 

Ek = k%Q.R, -N <k<N, (8) 

which are equidistant {AE^ = hflR), and again we con- 
sider that each spin interacts only with the field, individ- 
ually. 

III. QUANTUM TREATMENT OF 
SPIN-PHOTON COUPLING 



II. COMPLEX SUSCEPTIBILITY IN THE 
LINEAR RESPONSE THEOREM 

In the linear response theory^^'^^, the imaginary part 
of the complex susceptibility is given by 
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-/3uih 



/oo 
(M,(0)M,(t))oe-^"*dt, (5) 
-OO 



where = (• • •)o = Tr • • • e-'^^^/Tre" 

/3 = l/ksT, fee is the Boltzmann constant and T is 
system's temperature. This gives a coefficient of pro- 
portionality between the induced quantity (M^) and the 
field ho at ho = 0. The eigenvalues of Hzo are simply 
given by Ek = khojQ, k = —N, ■ ■ ■ ,N, each of which is 
nCu = N\/k\{N—k) \ times degenerate. There is only one 
energy difference AE = hoJo that has nonzero matrix el- 
ement of Mx- In this case, each spin interacts only with 
the field, individually. Thus, we have a single peak in 
x"(w) at w = Wg. If we include some interactions among 
spins, such as the dipole-dipole interactions, the degen- 
eracy of the energy levels of the N spin system would be 
resolved, and additional peaks in x" (w) are expected. 



Quantum dynamics of paramagnetic spins 
under an ac field: Rabi oscillation 



If we consider the dynamics of spins in the ESR Hamil- 
tonian "Hs, the total magnetization shows the so-called 



Now we study the case where the interaction between 
spins and photons is treated quantum mechanically. We 
start by reviewing the energy diagram of model (4). 



A. Case of single spin N = 1 

First, we consider the case of A'' = 1. We adopt the 

basis |n,(T) where n denotes the number of photons in 
the cavity, and a = —/+ shows the ground/excited state 
as an eigenvalue of cr^ . The matrix of H is separated into 
2x2 blocks for each pair {|n, — ), \n — 1, +)}, given by 



-hujr. 



%g\/n 



rigV^ ^ + (n ■ 



The eigenstates are given by 



l)hu 



E± = noj{n - 



(9) 



(10) 



where A = w — and 6 = A"^ + ing^. At resonance 
u) = LVo, the above reduces to 



E± — [n— ^)huJo ± hgy/n, 



(11) 



Note that all blocks have the same photon number -|- 
magnetization constant'', C = n — 1/2 = (n — 1) -|- 1/2. 
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B. Case of a spin ensemble AT > 1 



1. Case Un 







For TV > 1, the working basis becomes \n, {ai ■ ■ ■ (Tn}) 
where ffi ~ — /+ shows the ground/excited state of ith 
spin, as an eigenvalue of cr? . When photons are absorbed 
or emitted by the spin ensemble, the quantity C = n+M 
is conserved®^''', with M = m — N/2 the ensemble magne- 
tization, m the number of excited spins and n the number 
of remaining photons. The Hilbert space corresponding 
to all m values is nCq + . . .+n Cn = 2^ in length. Fur- 
ther division in independent sub-blocks can be done if 
one use a total spin representation^. The system Hamil- 
tonian can be written as: 

n = + hLjb% + hg{bS+ + b^S-), (12) 

with S^'~^'~ = 'l2i' ■ The Hamiltonian commutes 

with the total spin (S'^ f + {3^^ + {S^' f and can be 
further separated into smaller blocks, classified by the 
total spin S. When all spins are in the ground state, 
the total spin is maximum, S = N/2. From here on we 
discuss the case S = N/2, since further photon excita- 
tions and emissions will selectively couple states within 
this subspace only. The subspace can be further lim- 
ited if there are not enough photons to flip all the spins 
in the system: m ranges from to m.in{N , rimax) with 
nmax = n + m the number of photons for Sz = —N/2. 
Using a basis {|m),0 < m < min{N , nma.^)} and the 
relation {S = N/2): 

S+\S,M) = ^S{S +1)- M{M + l)!^, M+l) = 
= ^{m + l){N-m)\S,M + l), 

the diagonal and off-diagonal terms of Eq. (12) are: 

Um.m = (m - N/2)?iujQ + (rtmax - m)?iw 
y-m,m+l = ^V^^max - m^/ (m -|- 1)(A'' - m). 



(13) 
(14) 



At resonance, the diagonal term becomes Hm,m = 
{nmax — N/2)TujOo, and is independent on m. 



IV. EIGENVALUES FOR AN ENSEMBLE OF 
SPINS 

For a total spin S = N/2, one expects A''-|-l eigenstates 

of which N — Umax are diagonal spin states and rimax + 1 
are coupled spin-photon states (if rimax > N, all iV -|- 1 
states are coupled). 



A. Analytical expressions for rimax = — 3 



Analytical expressions for eigenvalues E't 
listed below for few simple cases. 



(^max) 



are 



. In the vacuum-field where no photon exist when all 
spins arc in the ground state, there is a unique state and 
its eigenvalue is given by 



(0) _ TiwqN 



2. Case 



(15) 



. In this case there are two coupled states: |n = 1, m = 
N) and \n = 0,m = N — 1). The Hamiltonian of this 

block is 



huji 



o{—Y)-^^ gvN \ \n—l,m = 0) 
VN htuo (-f + 1) ; |n = 0,m = 1) 



and the eigenenergies are given by (m = 0, 1): 



(16) 



Ei/^^/n = a;o (^-y + l) + ^ + (^m - y'A^ + ANg^ 

(17) 

The VRS is given by 

e[^^ - = ?iv'A2 + 4Ng-^ (18) 

and represents the rate at which the spin system coher- 
ently exchanges one photon with the radiation field. 



3. Case rimax = 2 

. The spin- photon states span over three states: \n = 
2,m = N), \n = l,m = N - 1), and \n = 0,m = N -2). 
At resonance A = 0, the Hamiltonian is given by: 



hu)o{2- 



N 







ngV2N 
hgV2N hg^2{N-l) 



ng^J2{N -I) 



(19) 



and the eigenvalues are given by (m = 0, 1, 2): 

E^^^ /^i = wo ^ + 2^ + (m - l)gVW^ (20) 

Here, the Rabi frequencies are degenerate 

£;(2) _ ^(2) ^ £;(2) _ £:(2) ^ hgVW^. (21) 
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4- Case Umax = 3 

. Following a similar approach as for n„ 
Hamiltonian at resonance is given by: 



2, the 







hcjo (3 - f ) + 

/ ngVm 

hg^/m ng^'i{N -I) 

ng^A{N - 1) hg^^N - 2) 

V ng^J^N - 2) 

(22) 

with eigenvalues: 



i;^V?i = ?iwo(-iV/2+3)±5^5(Ar - 1) ± v'(4iV - 5)2 + 87V, 

(23) 

where m = 0, 1, 2 and 3 counts the four possible values. 

B. General case iV, rimax > 1 

The size of the block representing the spin-photon 

states increases with rimax- For each photon made avail- 
able to the spin system, an additional spin will partici- 
pate in the; cooperative energy exchange and a new spin- 
photon state is generated, as indicated in Fig. 1(a). This 
effect is exemplified in Fig. 1, at resonance, for AT = 5 
where up to six states arc generated with increasing rimax- 
When rimax = 5, all spins are participating and the size 
of the matrix is bounded at 6 x 6. Although no new 
states arc generated for values of rimax > N , the + 1 
eigenvalues will adjust to indicate the "oversaturation" 
with photons (as shown with connected dots in Fig. 1). 
Analytically, this effect is shown by Eq. (10) for iV = 1 
or by replacing N with rimax in the off-diagonal terms 
leading to Eqs. (15), (17), (20), and (23) and thus in the 
corresponding Rabi splittings. 

For large values of nmax, all Rabi splittings 



E. 



(«„ax) 

m+1 



(24) 



are equal to « 2%^nniax (as shown by vertical arrows 
in the example of Fig. 1). However, this transition to 
equidistance is not a smooth process. In the following, 
we show that the spread in Rabi frequencies becomes 
maximal at rimax = and that a gradual transition 
towards an equidistant spectrum develops for rimax > N. 
Difference between consecutive eigenvalues are shown in 
Fig. 2. For iV = 10 [Fig. 2(a)] one observes a spread of 
Rabi splittings over few units of g for rimax = A'" = 10, 
followed by a collapse on a single-valued Rabi splitting 
2%Y^nmax [shown by continuous lines in Fig. 2(a) and 
(b)] for rimax larger than several tens. 
An exact diagonalization study for N up to 400 and 



three 



ralucs provides support for a general view of 



the process [Fig. 2(b)]. At rimax = N (black dots) devi- 
ations from the ^hg^rimax limit (black lines) are maxi- 
mal and the range of spread increases proportionally with 
VN. For rimax = N + 250 and N + 500 the distribution 




LU -20g 



FIG. 1. (a) For each additional photon, a new spin-photon 

state is generated, up to a number A'' + 1. (b) Spin-photon 
eigenstates for A*' = 5, at resonance, obtained from diagonal- 
ization of blocks limited by [miri(rimax, N) + 1], as sketched in 
the insert. For rimax ^ N, all A^-l-l states become equidistant, 
leading to equal Rabi splittings (shown by vertical arrows). 




FIG. 2. Rabi splittings calculated at resonance, as difference 
between consecutive eigenvalues, for (a) A^ = 10 and (b) N 
up to 400 and rimax = AT (dots), A" -|- 250 (squares), A^ -I- 500 
(triangles); the N = 1 limit for Rabi splitting, 2?i(/^nmax, 
is shown by continuous lines. At ramax = N the distribution 
width in Rabi splittings is maximal, and it decreases rapidly 
with increasing rimax. 
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width of Rabi splittings is highly reduced. A convergence 
toward 2%^nniax is observed especially at low values of 
N (or very large Uraax/N ratio), since in this limit one 
approaches the analytical case of = 1 [Eq. 11]. 

This convergence corresponds to the energy diagram of 
the Rabi oscillation in Eq. (8) in the ESR model, where 
each atom interacts with the field individually. The fact 
that fifl cx ^nmax in this limit, resides on the dependence 
ho oc ^nmax- If the number of photons is large enough, 
b, 6t ^ (&), ^ .JW;^ in Eq.(4). Thus, ha in Eqs. (2), 
(3) corresponds to y'rimax- 
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C. Photon transmission spectra: a quantum to 
classical transition 

To probe the spin-photon states one could use a 
low power beam and analyze the transmitted signal (an- 
other option, demonstrated experimentally in the last 
section, is to study the Fourier transform of the coherent 
emission of an excited spin-cavity system) . A low power 
probe could be used to scan the frequency response of 
the cavity, after the introduction of rimax photons of fre- 
quency oj (= Wo at resonance). 

A variable frequency beam, probing after the introduc- 
tion of rimax photons, will see the Rabi splittings given 
by Eq. (24). Therefore, in such experiment one should 
be able to detect a transition in peak distribution from 
a large number of fi/f values to a single-valued trans- 
mission peak. By design, such photon-driven transition 
cannot be observed &i N = 1 but it would require at 
least several non-interacting spins coupled to rimax ^ N 
photons. 

The passage of a system from Timax to rimax + 1 implies 
an excitation from one group of eigenstates to the next 
one, as sketched in Fig. 1(a). Due to the large number 
of eigenstates involved, the distribution of energy differ- 
ences 

?^AL"r^ = E^r-^'^ - ^i"--^ (25) 

can be broad. As shown in Fig. 3(a) for N = 10, the dis- 
tribution width is increasing as ^rimax- However, only in 
the vicinity of rimax = N are the values highly dispersed, 
as shown by a dense cloud of points at rimax « A^ in 
Fig. 3(a). For large values of rimax, and in full agreement 
with the study of Fig. 2, the A^,'J^'''' values are equally 
spaced by a Rabi splitting ^Ug^ rimax- 

A different probing method would consist in using the 
rimax photons to actually probe the energy levels. An ex- 
perimental demonstration for one superconducting qubit 
electrically coupled to rimax = 0...5 photons has been 
recently performed^® (see also Ref. 24 for a multi-atom 
experiment). In such a case, multi-photon transmission 
experiments can probe the energy difference between the 
and Eg [of Eq. (15)] 

huj = {E^-'^- - E°)/nmax. (26) 



s 

1 


2g- 


max 

single peak 

(classical) * * 




* Rabi splittings* 

(quantum) 






• 










• 




oo 
LU 


og- 




• 


• 


Je 


-2g- 


N=io '■■"i'-.. 

1 1 1 1 1 1 1 1 1 1 


• 

• 

• • • 

1 1 1 1 1 


(b) 

• • 

1 1 1 1 1 1 1 1 1 



0.0 0.2 0.4 ,0.6 0.8 



FIG. 3. (a) Representation of all possible excitations from 
rimax to rimax + 1, Calculated at resonance, for N = 10. 
The corresponding energy splittings are closely packed in the 
transition region rimax ~ N and became equally spaced by 
2?ig^nmax for rimax S> N. (b) Multiphoton, nonlinear reso- 
nances as a function of 1/-^ rimax. The linear dependence at 
high rimax indicates the emergence of an equally spaced Rabi 
spectra, and the convergence to a single peak indicates the 
gradual passage to a classical ESR, single-peaked, resonance 
condition. 

Above the transition threshold rimax ~ N, the photon fre- 
quency thus defined becomes gradually insensitive to the 
distribution in Rabi splittings due to the l/rimax factor. 
A level diagram with Rabi splittings uniformly spaced 
by 2%^nmax will generate transmission peaks spaced 
by 2hg / y/Hmax- The linear dependence on 1 / ^/rimax for 
rimax ^ N is shown in Fig. 3(b), calculated at resonance 
conditions for N = 10. The gradual passage from the 
quantum case to the classical ESR condition (single peak 
at w = Wo) is visible with the increase of the electromag- 
netic intensity. 

V. EXPERIMENTAL STUDY OF STRONG 
COUPLING IN SPIN SYSTEMS 

We demonstrate the strong coupling between an en- 
semble of spins s = 1/2 and photons using a sample 
of the well-known ESR standard materiaP^, dipheriyl- 
picri-hydrazyl (DPPH). For an optimized cavity cou- 
pling, sample positioning and size, we have been able 
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FIG. 4. (Color online) (a) Fourier transform of cavity ringing, 
measured at room temperature for several magnetic field val- 
ues around the resonance condition. Frequencies are relative 
to the sample-loaded cavity resonance, in absence of applied 
field. The observed two peaks are indicated by vertical marks, 
(b) Peaks position as a function of detuning field noSH^- The 
dashed lines indicate the classical level diagram, whereas the 
continuous line fit [Eq. 27] shows a measured Rabi splitting 
of 10.9 MHz. 



to induce a sizeable Rabi splitting observed for a series 
of Zeeman spin splittings luq. 

A l-/Lts-long microwave pulse pumps a large number 
of photons into a cylindrical cavity operated in mode 
TEgii. When the microwave is switched off, the cavity 
is coherently emitting photons corresponding to its own 
eigenmodes (phenomenon known as cavity ringing) . This 
ringing is detected by a homemade heterodyne analyzer. 
The experiment is performed at room temperature. To 
ease the distinction between the photons of the pump 
pulse and cavity's own emitted photons after pumping, 
the pump is detuned by 50 MHz from cavity's resonance 
uj/2tt — 9.624 GHz. The Fourier transform of the coher- 
ent oscillations is shown in Fig. 4(a), with the frequency 
axis shifted by uj/2tt for clarity. The sample-loaded, no 
field (wo = 0) oscillation shows only the cavity signature, 
located at w/27r. 

For an applied magnetic field fioHz = hoj/gsfis, with 
5s = 2, the spin system is in resonance with the cavity 
{ujQ — uj). The eight traces of Fig. 4(a) are for fields 
within ±0.4 mT of the resonance condition. The Fourier 
transform of the coupled spin-photon oscillations show 
peaks, indicated by vertical marks in Fig. 4(a). The 
peaks location, relative to cavity's resonance uj/2t:, is 
shown in Fig. 4(b), as a function of the field detuning 

In a classical ESR experiment, for instance using the 
DPPH as a field calibration standard, one expects energy 
levels that follow the dashed lines in Fig. 4(b). In par- 
ticular, the cavity peak is visible and changes abruptly 
in size but not location, when the resonance condition is 
met (at the intersection of the dashed lines). 



In our experimental conditions, due to the spin-photon 
strong coupling, one observes two peaks separated by the 
characteristic Rabi splitting. The size of the splitting can 
be attributed to the interaction between a group of N 
spins and a single photon, although the cavity contains 
a large number of nmax photons. A possible explanation 
resides on the existence, in standard ESR experiments, 
of the so-called "spin packets" characterized by a coher- 
ence time T2 and grouping an average number of spins N. 
The T2 and cavity decay times can be estimated from the 
peak widths of « 6 MHz and 2.7 MHz at uj = ujq (res- 
onance) and Wo = (no field) conditions, respectively: 
T2 = 170 ns and Tcav = 370 ns. 

Consequently, the cavity photon depletion during emis- 
sion, marks the transition between the Eq^I and E^^^^ lev- 
els of Eqs. (17), and (15). Cavity ringing shows coherent 
oscillations with the frequencies a;eo,i/27r, given by: 

uJeOA~^ = -^±l^A^ + ^% (27) 

on which the continuous lines of Fig. 4(b) are based. 
The fit procedure leads to a Rabi splitting of n/j/27r = 
10.9 MHz. We note that a detailed knowledge of the g 
and N parameters would require an on-chip type of ex- 
periment, at low temperatures, to maximize the "spin 
packet" size and to ensure a precise knowledge of spin 
position^^. In a single photon picture, and knowing the 
cavity volume (^ 50 cm'^), we estimate the number of 
spins contributing to flu to be on the order of 10^". Our 
data demonstrate the spin-photon strong coupling and 
also provide the mainframe to study experimentally the 
transition shown in Fig. 3(b). The role of other factors 
in the spread of Rabi splittings, such as: dipolar or hy- 
perfine interactions, local anisotropic crystal fields, and 
the size of the sample vs. size of a cavity (mode) can 
be studied as well. At the same time, the possibility 
to entangle spins and photons inside a cavity comes in 
strong support for the implementation of quantum com- 
puting algorithms by using magnetic spins and on-chip 
quantum electro-dynamics methods. 

In conclusion, we present a study of cooperative spin- 
photon interaction leading to a transition between a 
quantum-type Rabi spectra to a classical ESR spectra. 
The transition requires a number of spins iV > 1, and the 
possibility to gradually increase the number of photons. 
When the later is close to N, a dense Rabi spectra is nu- 
merically observed. For a large number of photons, the 
spectra gradually become equidistant and the resonance 
peaks converge toward the classical single-peaked, ESR 
resonance. The theoretical model is complimented by 
an experimental demonstration of the spin-photon strong 
coupling regime, which uses the B-field component of an 
electromagnetic mode in a cylindrical cavity. 
Note added. Recently, two similar experimental results 
have appeared^^'^^. 
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